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Abstract 

We calculate single pion production by neutrinos in the resonance region. 
We consider both charged and neutral current reactions on free protons and 
neutrons. We present differential and total cross sections which can be com- 
pared with experiments. Then we use these results to calculate the spectra 
of the emerging pions including the Pauli suppression factor and rescattering 
corrections for reactions in heavy nuclei. Our results will be useful for study- 
ing single pion production and for investigating neutrino oscillations in future 
experiments. 

I. INTRODUCTION 

There is strong evidence for the mixing of muon neutrinos with another state, being 
either tau or sterile neutrinos. The evidence comes from atmospheric neutrino experiments 
which observe a decrease of muon neutrinos in charged current reactions, but no decrease of 
the corresponding electron neutrino interactions 

In order to obtain better insight into the oscillation which takes place and in order to 
eliminate dependence on the flux there are proposals and experiments being planned and 
constructed, which look at the neutral current interactions. These are reactions which will 
use neutrinos of an average energy of 1 GeV producing the resonances between 1.0 and 
1.6 GeV/c 2 . One proposal considers the production of pions directly by the atmospheric 
neutrinos @ and the detection of 7r°'s with the help of two ring events. This method is 
limited by the low flux of atmospheric neutrinos. 

More powerful are experiments which use neutrinos from an accelerator with two detec- 
tors; the first one nearby the accelerator and a second further away. The nearby detector 
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will be able to detect all pions and check the flux, as well as the cross-sections for these 
reactions. The detector with the long-baseline (300-400 km) will observe the charged and 
neutral current reactions. 

We classify the reactions. Quasi-elastic scattering v^+n — > fi~+p is well understood. The 
oscillation of muon neutrinos into other neutrinos, tau or sterile, will produce a reduction 
of muon events in the far away detector. 

In addition to the above reaction there are excitations of resonances and their subsequent 
decays 



Vv + n — > jj, +n + 7r + 



Vp+tl — > fl + p + 7T° . 

Furthermore, there are the neutral current reactions 

_o 



Up, + p — > ^ + n + 7r+ 



z/ M + n — > Ufj, + n + 7T° 



1.1 



1.2 



1.3 



1.4 



1.5 



1.6 



1.7 



The theory for the production of these states is known for thirty years now and there are 
several calculations available. The charged current reactions have been studied extensively 
|||-f| and the production of the A ++ has been understood theoretically. It has also been 
measured experimentally with good agreement between theory and experiments. For the 
other charged current reactions there are few experimental measurements. For this reason 
the nearby detectors of the experiments should study the reactions using light and heavy 
nuclei as targets. 

Knowledge of the neutral current reactions is even more limited. Two of the latest 
calculations of charged and neutral current reactions for the production of single-pion cross 
sections differ by approximately 20% ||||. One should make all possible efforts now to 
reduce the overall uncertainty and measure the various channels experimentally. 

A second difficulty arises from the fact that the experimental targets are heavy materials 
so that the interactions take place on protons and neutrons bound in nuclei like for example 
sO 16 , 18^4?" 40 or 26-^e 56 . In the heavy nuclei, the produced pions rescatter before they exit 
from the nucleus and are subject to two phenomena: (1) the cross-sections are reduced by 
the Pauli exclusion principle, when the energy of the recoiling nucleon is low and can not 
occupy a filled level of nucleons, and (2) the pion charge exchange due to rescattering. These 
phenomena are known and have been the subject of extensive studies PJTOf. 
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The expectations of the experiments are the following. For all charge current reactions 
we expect a reduction of the observed rates in the far away detector because some of the 
muon neutrinos oscillated into another state. For the neutral current reactions we expect no 
reduction in rate if the oscillation is to tau neutrino because all neutrinos contribute equally 
to neutral current reactions. We expect a reduction if the oscillation is to sterile neutrinos. 

Now, since the reduction is expected to be approximately 40%, it is important to un- 
derstand all possible corrections. An important requirement is that the nearby and the far 
away detector use the same nuclei as targets. In case this is not possible, then corrections 
will have to be applied. 

Because of the importance of the experiments and the opportunities they present for 
establishing the charged and neutral current reactions, we have undertaken the task of 
calculating the cross-section on free protons and neutrons. This way we produce differential, 
as well as integrated cross-sections. Then we use the obtained results to calculate the 
corrections which are present in the nuclei. 

This paper is organized as follow: Section 2 is devoted to the calculation of the differ- 
ential and total cross-sections for single-pion production in neutrino-nucleon interactions. 
In Section 3 we discuss the nuclear effects involved in this process and calculate the energy 
spectra for charged and neutral pions for a few different materials typically used as targets 
in the experiments. Finally, results and conclusions are presented in Section 4. 



In this section we present the main equations and the form factors used to evaluate 
the differential and total cross sections for single-pion production in neutrino-nucleon in- 
teractions. For neutrino energies of a few GeV the single-pion production proceeds mainly 
through the excitation of the lower resonances. The main contribution to the cross section 
comes from the production and the subsequent decay of the A(1232)P33 resonance. Nev- 
ertheless, some of the channels receive a non negligible contribution from the isospin 1/2 
resonances as, for example, the N(1440)Pn and the N(1535)Sn resonances. 

The channels under investigation in this paper are the three charged current and the 
four neutral current channels listed in the introduction in Eq. 1.1 to 1.7. 

Using Clebsh-Gordan coefficient, it is easy to verify that the amplitudes for these seven 
channels are given by the following equations: 



II. SINGLE-PION PRODUCTION 



A. General Formalism 



A( f i-+p + n + )=A« 



(2.1) 



A(fT +n + n + ) = 1/3A£ + 2V2/3A? 



(2.2) 



A(fi~ + P + *"°) = -V2/SA™ + 2/3A c ; 



cc 



(2.3) 



A(u„ + p + 7r°) = V2/3A™ + \j3A\ c + 1/3A 







(2.4) 
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A{up + n + 7T+) = -1/3A" C + \/2/3^ c + y/2/ZA\ (2.5) 
A(z/„ + n + tt°) = y/2/3A? + 1/3A" C - 1/3A° (2.6) 

A{v„ + p + 7T-) = 1/3A" C - \/2/3Af + V2/?,A\ (2.7) 

where, in this paper, A^ c ' nc corresponds to the amplitude for the production of the P33 
resonance, Al c ' nc is the sum of the amplitudes for the production of the Pn and Sn resonances 
and A\ is the sum of the isoscalar contributions of the Pn and Sn resonances to the cross 
section. 

Notice that, as suggested in Ref. |§], the neutral current amplitudes A% c , A™ c and A\ 
can be derived from the corresponding charged current amplitudes A^ c and A^ by simply 
rescaling the vector and axial form factors. In the case of the A™ and A™ c amplitudes the 
vector and axial charged current form factors need to be multiplied by 1 — 2sin 2 8w and 
by 1 respectively, where 9w is the weak mixing angle. For the A\ instead the vector and 
axial charged current form factors need to be multiplied respectively by —2/3sin 2 8 w and 0. 
Furthermore, we want to point out that, since the A\ amplitude turned out to be very small 
compared to A^ c and A™ c , we neglected the isoscalar contribution in our evaluation of the 
cross sections. 



1. A(1232)P : 



As it was mentioned in the introduction, the theory for the production of the A(1232)P33 
is well known and understood, and several independent calculations have already been pub- 
lished, showing good agreement with the experimental results Therefore, rather than 
developing our own formalism for this process, we decided to follow the article of Schreiner 
and von Hippel in Ref. 0. 

The fully differential cross section da/dQ 2 dWdQ 7r from Ref. has been integrated over 
the polar angle 0^ and converted into the triple-differential cross section da/dQ 2 dWdE n 
by using the fact that E l £ b = ^E^ M + f3^\j% M \cos6 n , where the superscripts lab and CM 
denote quantities measured in the laboratory and in the center of mass frame respectively. 
The total cross section has been obtained by integrating over the allowed range of values for 
Q 2 and E n , with a cut on the invariant mass W range at 1.6 GeV. 

The axial and vector form factors used in this calculation are the ones given by Alvarez- 
Ruso et al. in Eq. 12, 13 and 18 of Ref. fIT |. Notice that, since these form factors have been 



derived from photo- and electro-production experiments in which a A + or a A was produced, 
in order to obtain the correct cross section for the A ++ production, all the form factors need 
to be multiplied by V3 due to the fact that < A++ | V Q | p >= < A+ | V° m | p >. 



2. N(1440)Pu and N(1535)S U 



As shown in Ref. [Q, the triple-differential cross section da/dQ 2 dWdE 7T for the produc- 
tion of the Pn and Sn resonances is given by the following equation: 
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da W GlQ 2 1 // i/\ lirlS If i>\ 2 , lr ,. . Iw 

//l MMJ 2 + -(1- — ] |M r | 2 + -|M,| 



dQ 2 dWdB w 4Mtv (27r) 4 i/ 2 y/v 2 + Q 2 \\ E V P 2 V 



(2.8) 



where M/y is the nucleon mass, Gp is the Fermi constant and v is the difference between the 
energies of the incoming and the outgoing lepton. Remembering that nearby the resonance 
only the s-channel is essential, the three matrix elements Mi with i — r,l,s can be defined 
as follow: 

Mi = -f R u{pi) lb (jt> + i + M R ){g v f - g A fas)u(p)f{W) (2.9) 

where /r is the coupling constant of the pion to the nucleon and the resonance, Mr is the 
resonance mass, p and pi are the four-momenta of the initial and final nucleon respectively, 
q is the four-momentum transferred from the leptons to the hadrons and gv,A are the vector 
and axial form factors. The values for are given in Eq. B7 and BIO in Appendix B of 
Ref. 1. 

The Breit-Wigner factor f{W) in Eq.2.9 can be written in the following way: 



where T R (W) = T R q w (W)/q 7T (M R ) with q n (W) = yj \W 2 - M% - M 2 ) 2 - AM 2 N Ml/2W and 
T° R the width of the resonance. 

Finally, the polarization vectors e l with i = r,l,s used in Eq.2.9 are defined as: 

e s = ^v* + Q*,0Av) , e^ = -J=(0, 1,^,0). (2.11) 

For the form factors we use the expressions obtained by Fogli and Nardulli in Ref. [Q. 

As in the case of the P33 resonance, the total cross section has been obtained by inte- 
grating over the allowed range of values for Q 2 and E n , with a cut on the invariant mass W 
range at 1.6 GeV. 



B. Results 

In this section we present our results for the total cross section for the seven channels 
under examination and, where possible, we compare these results with experimental data. 

In Fig.l the total cross section for the + p — > \T + p + 7r + process has been plotted 
versus the incoming neutrino energy. The data points have been taken from Ref. [|13|] (solid 



circles) and from Ref. (empty circles). As it can be seen, the agreement between the 
theoretical curve and the experimental data is quite good. 

Fig. 2 and Fig. 3 display the total cross sections for the z/ M + n — > + p + 7r° and 
Vp, + n — > + n + 7r + processes respectively, again plotted versus the incoming neutrino 
energy. In this case the data points have been taken from Ref. |]TB[ (solid circles), from Ref. 
fl4| (empty circles) and from Ref. [|15j (crosses). Also in this case the agreement between 
the theoretical results and the data points is reasonably good. 
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The difference between the theoretical and the experimental results can be partially ex- 
plained by taking into account the fact that, while the theoretical curves have been estimated 
imposing a cut on the invariant mass W at 1.6 GeV, the experimental points have been ob- 
tained without any cut. Notice also that we didn't include any non-resonant background in 
our evaluation of the cross sections. 

In the case of the neutral current interactions, the experimental results are presented 
in the form of ratios between each of the neutral current channels and one of the charged 
current channels. For these reason, Fig.4, Fig.5, Fig.6 and Fig. 7, which display respectively 
the total cross sections of the +p — > u^+p + ir , +p — > + n + n + , v^ + n — > u^ + n + n 
and Vp + n — » u^ + p + ir^ processes plotted versus the incoming energy, have no data points. 
Nevertheless, we compared our results with the experimental ratios from Ref. fT6|-|T9| and 
found that there is a reasonable agreement, even if, in some cases, the ratios measured by 
the different experiments differ a lot one from the other. 



III. NUCLEAR EFFECTS 

In Section 2 we discussed the reaction v + N — > I + N' + 7r ± '°, where N is a free nucleon 
(proton or neutron). In order to investigate the nuclear effects taking place in the exper- 
imental targets (for example, s^ 16 , i&Ar i0 or 2 6-^e 56 ), we need to study the modifications 
necessary for the reaction v + T — > / + T' + 7r ±,ci , where T is the nuclear target and T' is an 
unobserved final nuclear state. 

We visualize the reaction as a two step process with the neutrino interacting with individ- 
ual nucleons producing single pions and excited nuclei. The production process is influenced 
by the Pauli principle and the Fermi-motion of individual nucleons. The subsequent jour- 
ney of the pions is a "random-walk" of multiple scattering until the pion escape from the 
nucleus. In the multiple scattering the pions can exchange their charge. These phenomena 
have been studied in the past || and we shall adopt the formalism in order to calculate the 
energy spectra of the pions. We give enough details so that the reader has an overview of 
the model but for more details he or she should consult the original article. 



A. Charge Density Distributions 

Following Ref. ||, we treat the target as a collection of independent nucleons which are 
distributed in space accordingly to a density profile determined through electron-nucleus 
scattering experiments. For the charge density profile of sO 16 we adopt the harmonic oscil- 
lator model in which the density is given by: 

p(r) = p(0) exp(~r 2 /R 2 ) [l + C^ + C, (^) | (3.1) 



where R = a/K with K = y3(2 + 5C)/2(2 + 3C) and a the root mean square radius. 

For isAr 40 and 2eFe 56 we use the two parameters Fermi model and write the charge 
density in the following way: 

p(r) = p(0) (1 + exp ((r - CyCi))" 1 . (3.2) 
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The different parameters used in Eq. [3. 1,3.2] are given in Ref. [2(| and are summarized 
in Table 1. 



B. Charge Exchange 



In the first step the neutrinos interact with the bound protons and neutrons with these 
reactions allowed provided that the energy of the recoiling nucleon is above the Fermi sea. 
This brings a correction factor calculated in ||. In the second step the pions rescatter several 
times until they reach the surface of the nucleus and escape. It is important to notice that 
the process taking place during the rescattering depend only on the properties of the target 
nucleus and are independent of the leptons involved in the first step. The differential cross 
sections for leptonic pion production on nuclear and on free nucleon targets are related to 
each other through the so called charge exchange matrix M in the following way: 



M 



where 



/ da{ z T A - +) \ 

dQ 2 dWdE w 
da( z T A ;0) 

dQ 2 dWdE w 
da(zT A ; -) 

V dQ 2 dWdE n ) 



da(N T ; ±0) _ da(p; ±0) 
dQ 2 dWdE„ dQ 2 dW^L^ 



( da(N T -+) \ 

dQ 2 diyd^ 
da(iV T ; 0) 

dQ 2 dWdE v 
da(N T ; -) 

V dQ 2 dWdE w J 



(3.3) 



+ {A-Z) 



da(n; ±0) 
dQ 2 dWdE v 



(3.4) 



Its eigenvalues define beams of pions of specific charge combination, which offer a scat- 
tering they produce a beam, which is decreased by the appropriate eigenvalue. The complete 
scattering phenomenon is characterized by three transition probabilities f(A) corresponding 
to the three eigenvalues. They describe the probabilities of beams with eigenvalues A = 1, 
| and | to survive and exit the nucleus. The "random-walk" of the pions is a stochastic 
process and several solutions were found in Ref. 0. We collect here the main formulas for 
the calculation. When we transpose the final pion to the initial state, we obtain the system 
i\i + nf whose total isospin can be 0,1 and 2. We use this property to parameterize the 
charge exchange matrix in terms of three functions A, c and d 



(l-c-d 



M = A 



\ 



d 

c 




The overall factor A is given by 



A = g(W,Q 2 )f(l) 



(3.5) 



(3.6) 



with g(W,Q 2 ) being the Pauli suppression factor and f(l) the transmission coefficient for 
the state with eigenvalue 1. The others two functions are 



1 1 /(5/6) 1 /(1/2) 
3 2 /(l) 6 /(l) 



(3.7) 
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As mentioned already f(A) contains the dynamics of the multiple scattering for the A 
eigenvalues. Both the Pauli factor and several models f(A) were presented in Ref. @. It was 
also shown that approximating the multiple-scattering with scatterings in the forward and 
backward directions provides a very accurate approximation. Thus what is important is the 
effective profile of the nucleus that the pions sea. This allows one to write 

fm f °°dbbL(b)f(\,L(b)) 

where b is the impact parameter and the effective length L(b) is given by: 

Lib) = JLj^dzpy/(z* + V) . (3.10) 
In the case of sO w , the effective length L(b) is given by: 

Lib) = RV^exp-b 2 /R 2 (l +c(± + ^J (3.11) 
while for i$Ar and 26^e 56 L(b) is written as: 



The appropriate expression for the function /(A, L(b)) has been derived in Ref. || both for 
the case of a one- dimensional multiple scattering problem and for the case of a spherical 
one. The two solutions have then been compared showing excellent agreement over the entire 
range of parameters. Therefore, in this paper, we adopted for f(X,L(b)) the approximate 
expression obtained in Ref. for the one-dimensional problem. 



C. Averaging Approximation 

It is important to notice that, while the Pauli production factor depends on both W 
and Q 2 , the function /(A) depends only on W and this dependence is very weak. Therefore, 
as it has been already verified in Ref. 0, it is reasonable to average the charge exchange 
parameters over the leading IV-dependence by defining an averaged function /(A) in the 
following way: 

JdWg(W0-V 3 ,3(W0/(A) 
/(A) " fdWq(W)-ia 3 , 3 (W) 

where 0^3 (W) is the pion-nucleon scattering cross section and q(W) is the pion momentum. 
For the definitions of cr 3i3 (W) and qiW) see Appendix C of Ref. 0. 



S 



D. Results 



Using the model outlined in the previous subsections, we evaluated the nuclear correc- 
tions for leptonic pion production on three different nuclei: oxygen, argon and iron. 

The values obtained for /(A) and for g(W, Q 2 ) for the three different targets under 
consideration are listed in Table II and Table III, respectively. Notice that, as Table III 
shows, the reduction of the cross section due to Pauli exclusion principle is larger for smaller 
values of Q 2 and it does not depend on the material. 

The charge exchange matrices M for oxygen, argon and iron are given by: 

/ 0.782 0.161 0.057 \ 
M( 8 16 ) = A 0.161 0.677 0.161 (3.14) 
V 0.057 0.161 0.782 / 



M( 18 Ar w ) = A 



I 0.731 0.187 0.082 
0.187 0.625 0.187 
\ 0.082 0.187 0.731 



(3.15) 



M( 26 Fe 56 ) = A 



/ 0.718 0.194 0.088 
0.194 0.612 0.194 
V 0.088 0.194 0.718 



(3.16) 



The Pauli factor and the charge exchange matrix M for oxygen have been compared 
with the corresponding quantities previously evaluated in Ref. P,|T0"1| and have found to be 
in good agreement with each other. Unfortunately, no comparison with previous calculation 
or experimental data is possible for argon and iron. 

The differential cross sections evaluated in Sect. 2 for free nucleon targets have been used 
here together with the charge exchange matrix M to obtain the differential cross sections 
for nuclear targets. These cross sections have been integrated over W and Q 2 keeping the 
neutrino energy fix at 1 GeV in order to obtain the pion energy spectra appearing in Fig. [8- 
16]. 

Fig. [8-10] display respectively the pion energy distributions for positive, neutral and neg- 
ative pions produced on oxygen targets. In each figure the solid line represent the initial 
distribution without any nuclear correction, the dashed line represents the same distribution 
after the application of the Pauli factor in the production, and the dotted line represents 
the final distribution after applying all the nuclear corrections discussed in the previous sub- 
sections. Similarly, Fig. [11-13] and Fig. [14-16] display the corresponding pion distributions 
produced on argon and iron targets, respectively. 

From these figures it is clear that, while the reduction of the cross section due to the Pauli 
production factor is the same for all the processes investigated in this paper, the nuclear 
corrections related to the pion charge exchange and pion absorption are larger for neutral 
pions than for the positive or negative ones. Furthermore, these corrections turn out to be 
larger for heavier nuclei. Finally, the magnitude of the nuclear corrections decreases with 
increasing pion energy. 
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IV. CONCLUSIONS 



At least three long-baseline neutrino experiments plan to study low energy neutrino 
reactions. Their main aim is the observation and better understanding of the neutrino- 
oscillations, but a necessary input is the understanding of these reactions in free protons 
and neutrons, as well as the modifications brought about when the nucleons are bound in 
relatively heavy nuclei. 

In order to work in a coherent framework we calculated the cross sections on free protons 
and neutrons. The theory for the production of the A(1232)P 33 resonance is well understood 
and our results for the total cross section agree with the experimental data. The same holds 
for other channels where 1=1/2 resonances also contribute. The comparisons appear in 
Figures 1-3, where it is evident that the accuracy of the measurements is subject to large 
improvements. Thus it is highly desirable that the new experiments use the nearby detector 
in order to measure the various cross sections. This refers to charged and neutral currents 
interactions on free protons and neutrons. The main uncertainties on this part of the paper 
are the functional form and parameters of the form factors and interference between 1=3/2 
and 1=1/2 resonances. We expect that the effects from these uncertainties are small. 

More important are changes which are brought about in the scattering of neutrinos in 
heavy nuclei. It is very likely that the far away detectors will use heavy materials as targets in 
order to enhance their counting rates. The heavy materials bring in corrections comparable 
to oscillations. In this article we used an old model for nuclear corrections and calculated 
the effects on the produced 7r ± '°. In Section 3 we reviewed the main features of the model 
so that the interested reader can reproduce the results. 

We decided that an interesting and important parameter in the experiments is the energy 
of the emerging pion. We calculated in Figures 8-16 the pion spectra as function of their 
energy. We found that the largest correction appears in the spectrum of the 7r°'s. The 
reduction of the signal for neutral pions with energies approximately 200 MeV is substantial: 
of the order of 40%. Processes with nuclear corrections as large as the ones found in this 
article require special attention. Several strategies suggest themselves. 

One is to use the same material for the front and the far away detector and study the 
spectra as a function of E n . Then compare the results from the two detectors and with quasi 
elastic scattering. In case that the experiments are forced to use different materials detailed 
calculations for the two materials will point to similarities and possible differences between 
the two targets. 
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TABLES 



TABLE I. Charge density distribution paprameters 



zT A 


a[fm] C[fm] C\[fm] R[fm] 


p(0)[fm- 3 } 


sO 16 


2.718 1.544 1.833 


0.141 


18 Ar i0 


3.393 3.530 0.542 4.830 


0.176 


26 Ee 56 


3.801 4.111 0.558 4.907 


0.163 




TABLE II. /(A) with A = 1, 1/2 and 5/6 for 8 O ie , i 8 ^r 40 and 


26 Ee 56 . 


zT A 


/(I) 7(1/2) 


7(5/6) 


8 o 16 


0.811 0.418 


0.587 




0.648 0.284 


0.420 


26 Ee 56 


0.625 0.261 


0.393 




TABLF, TTT Panli nrnHnrtinn fartnr n(W D 2 1 evaluated fnr W — 


l.2GeV. 


Q 2 




g(1.2GeV,Q 2 ) 


0.00 




0.785 


U.UO 




U.o 1 1 


0.10 




0.926 


0.15 




0.960 


0.20 




0.980 


0.25 




0.992 


0.30 




0.997 


0.35 




1.000 


0.40 




1.000 
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FIGURES 




FIG. 1. Total cross section for the + p — > \i + p + n + process plotted versus the incoming 
neutrino energy. The data points are from Ref.[13] and Ref.[14]. 
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FIG. 2. Total cross section for the + n — > \i + p + ir° process plotted versus the incoming 
neutrino energy. The data points are from Ref.[13], Ref.[14] and Ref.[15]. 
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FIG. 3. Total cross section for the Vn + n — > + n + 7r + process plotted versus the incoming 
neutrino energy. The data points are from Ref.[13], Ref.[14] and Ref.[15]. 
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FIG. 4. Total cross section for the + p — > + n + ir + process plotted versus the incoming 
neutrino energy. 
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FIG. 5. Total cross section for the + p — ► + p + 7r° process plotted versus the incoming 
neutrino energy. 
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FIG. 6. Total cross section for the + n —>■ + n + ir° process plotted versus the incoming 
neutrino energy. 
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FIG. 7. Total cross section for the f^ + n— >f (U +p + 7r process plotted versus the incoming 
neutrino energy. 
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FIG. 8. Pion energy distribution for positively charged pions produced on oxygen targets. 
The solid, dashed and dotted lines represent respectively the pion energy distribution without 
any nuclear correction, including only the Pauli production factor g and including all nuclear 
corrections. 
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FIG. 9. Pion energy distribution for neutrally charged pions produced on oxygen targets. The 
solid, dashed and dotted lines represent respectively the pion energy distribution without any nu- 
clear correction, including only the Pauli production factor g and including all nuclear corrections. 
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FIG. 10. Pion energy distribution for negatively charged pions produced on oxygen targets. 
The solid, dashed and dotted lines represent respectively the pion energy distribution without 
any nuclear correction, including only the Pauli production factor g and including all nuclear 
corrections. 
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FIG. 11. Pion energy distribution for positively charged pions produced on argon targets. 
The solid, dashed and dotted lines represent respectively the pion energy distribution without 
any nuclear correction, including only the Pauli production factor g and including all nuclear 
corrections. 
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FIG. 12. Pion energy distribution for neutrally charged pions produced on argon targets. The 
solid, dashed and dotted lines represent respectively the pion energy distribution without any nu- 
clear correction, including only the Pauli production factor g and including all nuclear corrections. 
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FIG. 13. Pion energy distribution for negatively charged pions produced on argon targets. 
The solid, dashed and dotted lines represent respectively the pion energy distribution without 
any nuclear correction, including only the Pauli production factor g and including all nuclear 
corrections. 
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FIG. 14. Pion energy distribution for positively charged pions produced on iron targets. The 
solid, dashed and dotted lines represent respectively the pion energy distribution without any nu- 
clear correction, including only the Pauli production factor g and including all nuclear corrections. 
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FIG. 15. Pion energy distribution for neutrally charged pions produced on iron targets. The 
solid, dashed and dotted lines represent respectively the pion energy distribution without any nu- 
clear correction, including only the Pauli production factor g and including all nuclear corrections. 
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FIG. 16. Pion energy distribution for negatively charged pions produced on iron targets. The 
solid, dashed and dotted lines represent respectively the pion energy distribution without any nu- 
clear correction, including only the Pauli production factor g and including all nuclear corrections. 
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